Plasmon electron-hole resonance in epitaxial graphene 
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Abstract 

The quasiparticle dynamics of the sheet plasmons in epitaxiahy grown graphene layers on 
SiC(OOOl) have been studied systematicahy as a function of temperature, intrinsic defects, in- 
fluence of multilayers and carrier density. The opening of the inter-band loss channel appears as a 
characteristic upward shift in the plasmon dispersion and a dip in the width of the loss peak, which 
is explained as a resonance effect in the formation of electron-hole pairs. Despite the existence of 
strong electronic correlations, the plasmon dispersion can be quantitatively described within the 
framework of a nearly free electron gas. 

PACS numbers: 81.05.ue, 79.20.Uv, 73.20.Mf 
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Graphene, as a two-dimensional (2d) lattice of sp^ -hybridized carbon atoms, proves to 
be a model system to study fundamental questions regarding electronic correlations, col- 
ective phenomena, many-body interactions, dynamical processes and their interrelations 
-|3|. Single layers can be synthesized by different techniques, e.g. exfolation of graphite or 
decomposition of hydrocarbons on transition metal surfaces [4,5]. Well orientated graphene 
layers on insulating substrates can be grown on SiC substrates by sublimation of Si. While 



the morphology, the interface, the e 
been studied extensively in the past 



ectronic structure, and even transport properties have 



recently 



I1J13 



10| , plasmons in graphene have attracted notice only 



13|. As we will show, for studies of this prototype of an excited collective state 



the prospects of plasmons of 2d films or even 
IigI . The comparatively fiat dispersion in these 



in 2d, graphene is indeed a perfect model system. 

Although plasmons are used in nano-scaled systems, e.g, for opto-electrical conversion 
and resonant energy transfer processes [m , 
of atomic wires are still at its beginning [15, 
systems provides ultrashort plasmon wavelengths. Thus, confinement on the nanometer 
scale becomes feasible, but strong damping of plasmon excitations and their sensitivity 
to defects on the atomic scale are challenges to cope with. Therefore, these studies require 
geometrically perfect and at the same time electronically flexible materials, such as graphene. 

On the other hand, the present knowledge on Id and 2d sheet plasmons is quite limited 
and, in particular, the influence of defects and the interaction with other quasi-particles is 
barely understood. It is clear that plasmon damping becomes dominant once the plasmon 
dispersion curve intersects the single particle band leading to Landau damping by single 



particle intra-band excitations (SPEintra) as seen, e.g., in Ag/Si(lll) 



17 



isj. Outside the 



Landau regime, this mechanism is still effective, but cannot be described without explicit 



considerations of correlations 



19l-|21| or momentum transfer processes [12|. Energetically low 



lying valence band states are considered only in terms of a polarizable background causing 
a quasi-linear dispersion at long wavelengths [l5| . However, couplings to other particles are 



neither included in random phase based descriptions 



20|, nor were signatures of resonant 



coupling seen so far in experiments. Only for hybrid systems, where plasmons interact with 



molecular orbitals or excitonic states 



22| . a resonant coupling was recently reported. 



In this Letter the sheet plasmon dispersion of graphene and the coupling to various 
decay channels is analyzed in detail by angle resolved high resolution electron energy loss 
spectroscopy (EELS). As we will show, once inter-band transitions are allowed by energy 



and momentum conservation {SPEmter), resonant coupling of this new decay channel leads 
to pronounced de-tuning of the plasmon frequency around Fermi energy and to an enhanced 
damping. The relevance of a similar type of quasiparticle dynamics in graphene, namely 
resonant coupling between plasmon and hole, called plasmaron, during photoemission was 
recently observed by Bostwick et.al. [2^. It corresponds to a complementary process to 
that studied here, and demonstrates the importance of resonant couplings for the plasmon 
decay in graphene. Thus, our findings broaden the standard picture for plasmon excitations 
and damping mechanisms. 

The growth of graphene and the measurements were performed under ultra high vacuum 
(UHV) condition. As substrate. Si-terminated 6H-SiC(0001) samples (n-doped, ~ lO^^cm"^ 
from SiCrystal AG) were used and graphene was grown by sublimation of Si while annealing 
the SiC substrate to approximately 1500 K. Further details can be found in Refs. jol and 
I2I . The plasmons were measured by using a combination of high resolution electron loss 



spectrometer (EELS) as electron source and detector with a low energy electron diffraction 



system (LEED) providing simultaneously high energy and momentum (A;||) resolution [24 1. 
Typical operating parameters were 25 meV energy resolution at a k\\ resolution of 1.3 x 

10-2 A-^ 

A sequence of EEL-spectra of a monolayer (ML) of graphene taken at different fcy values 



along the F-K direction is shown in Fig. [TJ 
attributed to the sheet plasmon in graphene 
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he broad loss peak can unambiguously be 



12| . A detailed analysis of the loss peaks (see 



Fig. [Hand Ref. [12|) reveals that the loss spectra can be decomposed into the main plasmon 
loss, two low energy phonon losses at 70 meV and 150 meV without significant dispersion, 
and small contributions from higher excitations, e.g, multipole plasmons. Already from 
Fig. H^) it is obvious that the plasmon losses reveal two main characteristics: the halfwidth 
of the loss, already large at low k\\, suddenly increases strongly around 0.1 A^^ (cf. with 
Fig. H]). Secondly, also the peak position shifts more strongly in this range of momentum. 
Both features are directly linked to the opening of a new decay channel by inter-band e-h 
excitations, as we will show. 

For this purpose, we first discuss the plasmon dispersion, shown in Fig. [21 revealing a 
dependence at small k\\. It rises almost linearly up to kp = 0.07±0.01 A"^, but then increases 
strongly again at higher k\\ values, very much in contrast to other known plasmon dispersions 



in 2d systems 



13, 



25| . This dip in the plasmon dispersion is an intrinsic property of pristine 
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FIG. 1. (color online) a) Loss spectra (semi-log scale) for 1 ML graphene grown on H-etched 
SiC(OOOl) as a function of momentum transfer, Primary electron energy 20 eV. Only the 
contribution of the 2d sheet plasmon to the loss spectrum is highlighted, b) Details of the fitting 
model, shown exemplarily for k = 0.05 A^^ , containing, apart from the main plasmon peak (solid 
line (orange)), a Drude tail (dashed (green)), phonon- and multiplasmon losses (dotted lines). For 



details of the fitting procedure see Ref. 
SiC(OOOl) (E=140eV). 



. c) LEED pattern of 1 ML of graphene grown on 



graphene, and does not depend on defect concentration, coupling to graphene multilayers 
or on temperature. Therefore, as seen from the coincidence of the dispersions measured at 
80 K and at 300 K (Fig. an increase of the carrier concentration by thermal activation 
can be ruled out as well as the coupling to high-frequency phonons. Although defects, 
e.g., atomic steps, influence sensitively the lifetime of the plasmons, particularly outside 
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FIG. 2. (color online) Plasmon dispersion for 1 ML graphene layer (average terrace length 150 nm) 



grown on H-etched samples, measured at 300 K (Q' 



free electron gas {N = 1 x 10 cm , 0.06 m, 



and80K(n). Fitted theoretical curves: nearly 



29( 1) and with dynamical screening 



20] (dashed- 



dotted). Dotted lines: Boundaries for intra- and inter-band SPEs derived from the graphene band 
structure. Inset: Independence of dispersion on step density (examples shown with terrace widths 
of lOnm and 150 nm) and on interaction to multilayers. 



the Landau damping regime 



12j, the characteristics of dispersion is again unaffected, as 



exemplarily shown in the inset of Fig. [2J Even in presence of multilayers the dip is visible 
and unshifted. Only for large k-values a small redshift due to the high polarizability of 
the graphene layer underneath is found, which is consistent with RefJlSl. The position of 
the dip in the dispersion does neither coincide with step periodicities nor with potential 



modulations caused by the superlattice of the bufferlayer underneath 



26|. Hence, so-called 
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zone boundary collective state effects [27|], as seen, e.g., in Al films [28!], are not crucial in 
graphene. Since the filling factors are comparably low, also intervalley scattering between 
Dirac cones can be excluded 2l| . The dip position, however, is very close to the Fermi wave 
vector at the given doping concentration {Ep ^ AOOmeV). 

Interestingly, existing theoretical models describe only sections of this dispersion. E.g., 
the nearly free electron gas model (NFEG) with inclusion of first order non-local field 



effects, as derived by Stern 29|], describes the experimental data well in the limit of large 
momentum transfer (k > 0.1 A"^) assuming an electron density of = 1 x and 
an effective mass of m* = 0.06 ± 0.01 rrie (solid line in Fig. [2]). These values are close to 
those deduced from photoemission [61] . Inclusion of dynamic screening for both electrons and 
holes (dashed-dotted line, shifts the dispersion to lower values and now yields a perfect 
match to experiment for small k\\ up to 0.09 A~^, but underestimates the plasmon energies 
in the SPEmter regime. Although e-h excitations are included in these RPA descriptions, 
the NFEG models do not contain effects like resonant energy transfer with corresponding k- 
dependent changes of lifetimes. As we will show below, only the inclusion of resonant energy 
transfer at Ep due to opening of the inter-band e-h excitations channel allows a quantitative 
fit of the measured dispersion curve, now even within an extended NFEG model. 

In order to observe such resonance effects between excited states, a short, but still suf- 
ficiently high lifetime Tf,_h at resonance is necessary in order to discriminate both types as 
separate excitations. As sketched in Fig. [3)d) the finite value of Ep above the Dirac point 
in pristine graphene enables inter-band transitions at finite k-values, which are shown as 



vertical transitions between shifted conduction and valence band states [30|. Indeed, we 
can describe the whole dispersion curve in every detail (see Fig. [3ti) from the loss function 
(Jm(— i/e)), using simply the Drude model of the dielectric function e = 1— '^p/w(w+i7), and the 
NFEG model of Stern assuming a resonance by opening the SPEinter channel. This results 
in a momentum-dependent dissipation '~f{k) for the plasmons. The effective damping '~f{k) is 
composed of a non-resonant term jnon-res, which takes into account intra-band transitions 
and structural imperfections, and a resonant part 7r.es. Fig. [3b) shows the best fit result for 
7. While '^non-res iucrcascs linearly with momentum 13], 'jres has a pronounced maximum 
around kp, as expected. The asymmetry in the damping curve demonstrates further that 
intraband transitions at low energies contribute to the dynamic loss channel as well. This 
pronounced damping is clearly seen in the experiment (Fig. [3ll), although the increase of 




FIG. 3. (color online) a) Fits of the plasmon dispersion with the NFEG model including damping 
(/ia;5jerra,7(fc)) thick(green) line) and resonant plasmon electron-hole coupling {Ejj and El, coupling 
constant 2A = 100 meV). For A = or damping constant j{k) — > the upper line {hw stern^^y^o) 
is obtained, b) Band diagram of pristine epitaxial graphene at finite k-values demonstrating 
intra-band inter-band e-h pair formation at finite k\\ values, c) Effective damping ^{k) and the 
decomposition into a resonant (res.) and non-resonant (non-res.) part, d) The FWHM of the 
plasmon loss is in qualitative agreement with 7(/c), shown in c). For details, see text. 

FWHM is smaller than calculated. Finite k-resolution and the simultaneous excitations 
of, e.g, phonons, not considered here, may be possible reasons. Nevertheless, the resonant 
broadening is even seen in much broader spectra, e.g. in presence of F4-TCNQ molecules 
acting as defect scatterers. 

This purely phenomenological description of the dip can be rationalized by a simple two- 
oscillator model representing plasmon and e-h-pair formation by inter-band transitions put 
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on-top of the non-resonant decay mechanisms. It yields a description fully consistent with 
the first approach. In order to see this de-tuning effect in the dispersion, the lifetimes, r, of 
both oscillators need to be of same order at kp . If further the interaction strength 2 A ^ fi/r, 
polaritons are expected in this regime of strong coupling. The upper and lower polariton 



states are given by ^c/,l(A;) = ^/2{Epi{k) + Eq) ± ^h^{2Af + {Epi{k) - Eo)^ where Epi and 
Eq are the non-interacting plasmon mode and exciton energy, respectively 3l|, |32|. Here we 
set the sheet plasmon energy hwstem for Epi and the Fermi energy Ep for Eq, which is the 
lowest energy of interband e-h pair formation without additional momentum. 2A measures 
the interaction between the quasi-particles. The best fit yields an interaction energy of only 
100 meV, in qualitative agreement with Ref. jsl. Consistent with the strongly enhanced 
damping determined from the fits described above, the observation of separate Rabi split 
states is suppressed, and only the continuous transition from the lower {E^) to the upper 
branch (Eu) can be observed. In fact, also this model fits the data quantitatively with values 
of damping close to those determined above. Formally, this model of plasmon electron-hole 
interaction is closely related to interacting plasmons located in spatially separated graphene 
layers which shows similar effects in the dispersion. However, since smooth single 
graphene sheets were used here, we conclude that the coupling to e-h oscillators is the 
dominant process. 

Finally we show that the position of the dip in k-space is only influenced by the chemical 
potential rather than by structural parameters. For this purpose, we adsorbed a small 
amount of F4-TCNQ 33|, which reduces the carrier density and hence shifts Ep downwards. 
Fig.Hlshows the dispersion after adsorption of approximately 0.01 ML at 300 K. As expected, 
the dip shifts to lower /cy-values (~ 0.03 A^^). This shift is also seen in the resonance 
of the FWHM. Since the added molecules act as scattering centers, the loss structure is 
further broadened so that the exact determination of kp from the FWHM is more difficult. 
Remarkably, the average slope of the dispersion remains unchanged, although the electron 
density at Ep is reduced to 3 x lO^^cm"^ due to adsorption of F4-TCNQ, which may have 
two possible reasons: The effective mass may be reduced by approximately the same factor 
when reducing the filling of the bands above the Dirac point, leaving N/m approximately 
constant. Secondly, the coupling of the plasmon mode with the loss channels, as evident 
from the large FWHM, leads to an effective integration over the electron density around 
Ep with a width proportional to the measured FWHM of the plasmon mode. Both the 
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FIG. 4. (color online) Graphene plasmon dispersion before (Q) ^-nd after F4-TCNQ doping (□). 
Inset: corresponding changes of k-dependent FWHM. Upper inset: stereographic model of F4- 
TCNQ. Solid line: Fit with the NFEG model without resonant damping. 

dynamics of the single and the collective excitation modes are extremely short and on the 
same time scale, which results in mixing of collective and single particle excitations. 

In summary, the detailed analysis of plasmon excitations and their decay mechanisms, 
using graphene as a 2d model system, yields ultra-short dynamics with life times of the 
order of r pa 10~^^s. Both e-h pair excitations by intra- and inter-band transitions con- 
tribute significantly to the decay of plasmons. When the latter channel is opened at finite 
k\\, coupling between plasmon and e-h pair excitations is resonantly enhanced. As a conse- 
quence, characteristic modifications of the plasmon dispersion and effective integration over 
the density states at Ep responsible for plasmon formation were found, thus widening the 
standard picture for plasmon excitations. 
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